Atmospheric neutrino oscillation analysis with sub-leading effects in 
Super-Kamiokande I, II, and III 
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We present a search for non-zero #13 and deviations of sin 2 O23 from 0.5 in the oscillations of 
atmospheric neutrino data from Super-Kamiokande-I, -II, and -III. No distortions of the neutrino 
flux consistent with non-zero #13 are found and both neutrino mass hierarchy hypotheses are in 
agreement with the data. The data are best fit at Am 2 = 2.1 x 10 -3 eV 2 , sin 2 613 = 0.0, and 
sin 2 #23 = 0.5. In the normal (inverted) hierarchy 613 and Am 2 are constrained at the one-dimensional 
90% C.L. to sin 2 6>i 3 < 0.04(0.09) and 1.9(1.7) x 10~ 3 < Am 2 < 2.6(2.7) x 10~ 3 eV 2 . The atmospheric 
mixing angle is within 0.407 < sin 2 6> 23 < 0.583 at 90% C.L. 
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PACS numbers: 14.60.Pq, 96.50.S- 



INTRODUCTION 



II. SUB-DOMINANT EFFECTS IN 
ATMOSPHERIC NEUTRINO OSCILLATIONS 



Despite experimental measurements of solar [THE], re- 
actor [9 , atmospheric [10] [TT , and accelerator [12] [13] 
neutrinos constraining their flavor oscillations, the na- 
ture of the neutrino mass hierarchy and whether or not 
#13 is zero remain open questions. The latter is the last 
unknown mixing angle and is currently the subject of a 
research program including beam and reactor-based ex- 
periments [14HT9] . At present, experiments have placed 
upper limits on the value of #13 [TO] I20H22] with the 
most stringent limit set by the Chooz [23] experiment. 
However, a non-zero value may manifest itself and be 
observable in the event rate of multi-GeV electron neu- 
trinos passing through the Earth, and to a lesser ex- 
tent, in similarly energetic upward-going muon samples. 
Though atmospheric neutrino data are well fit to pure 
v \i ^ v r oscillations with "maximal atmospheric mix- 
ing" [11 (#23 = tt/4), O v e transitions driven by 
solar oscillation parameters appear at sub-GeV energies 
when the atmospheric mixing deviates from this value. 
The questions of whether or not #23 is exactly 7r/4, the 
nature of #13, and the sign of the neutrino mass hierarchy 
all contribute to an eight-fold degeneracy [24] of oscilla- 
tion parameter solutions when considering CP-violation 
in neutrinos. For future experimental searches of CP- 
violation, answers to these questions are essential. 

In this paper two analyses are presented searching for 
evidence of sub-leading (second order) oscillation effects 
which address these questions and appear as changes in 
the v e and fluxes of the atmospheric neutrino sam- 
ples at Super-Kamiokande (Super-K, SK). The first is 
an improved extension of a three flavor oscillation anal- 
ysis using the first phase of the experiment (SK-I) [TO] . 
An updated analysis using the first, second (SK-II), and 
third (SK-III) phases is presented here. The data are 
then used in the second analysis to test whether #23 de- 
viates from 7r/4. 

The paper is organized as follows. In section [IT] we 
describe the oscillation framework used in the analyses. 



Section III discusses the data sample including additional 
sample selections designed to improve the sensitivity of 
each analysis. The methods and results of both are then 



presented in section IV and concluding remarks are found 
in section [V] 
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Neutrino oscillations in three flavors are described by 
six parameters: two mass squared differences, Am^, 
Am 2 3 , where Arn-^ = rrij— mf , a CP violating parameter 
(5 C p, and three mixing angles 9ij,(i < j). Each mixing 
angle parameterizes a rotation, [7^-, between mass states 
inside of the three-dimensional oscillation space. The cor- 
respondence between neutrino mass eigenstates and their 
flavor eigenstates is then: 



(i) 



1 

U = I c 2 3 s 2 3 

-s 2 3 c 2 3 



where U is the 3 x 3 PMNS matrix [25] [26] defined by 

^23^13^12, 

C13 s 13 e~ l6c 

1 

-s 13 e i6c p ci 3 

C12 S12 \ 
x I -5i2 C12 . (2) 
1/ 

Non-zero mixing angles and non-degenerate mass eigen- 
values give rise to standard neutrino oscillations. Obser- 
vations of solar and reactor neutrinos are well-described 
by oscillations governed by the "1-2" (solar) parameters 
while those of atmospheric and accelerator neutrinos are 
described by the "2-3" (atmospheric) parameters. These 
measurements have established two oscillation frequen- 
cies which differ by a factor of ~ 30. The third set of pa- 
rameters has been probed by Chooz, a reactor neutrino 
disappearance experiment sensitive to oscillations at the 
atmospheric Am 2 , which placed a limit on mixing in this 
channel at sin 2 6>i 3 < 0.04 for Am 2 - 2.0 x 10" 3 eV 2 at 
90% confidence 02]. 

For the purposes of studying sub-dominant oscillations 
in atmospheric neutrinos, it is useful to consider oscilla- 
tion probabilities in two domains: (i) #13 ~ such that 
U13 ~ I, and (ii) #13 > 0, but oscillations driven by the 
solar parameters are negligible. The observable effects 
each domain has on the atmospheric neutrino sample can 
similarly be divided into two energy regimes motivating 
two separate analyses with distinct foci: each analysis 
has been tailored to its regime of interest. 

In the case of #13 ~ 0, the neutrino oscillation proba- 
bilities in constant density matter may be written [28] 

P(v e = COS 2 2 3Pex (3) 

P(v» = 1 - cos 4 (9 2 3^ (4) 

-sin 2 2(9 23 (l - \A - Pexcoscj)) 
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FIG. 1: (color online). The left side of the figure shows the calculated v e transition probability P ex for atmospheric neutrinos 
with an energy E u and neutrino zenith angle, cosG^, using Ara? 2 = 7.7 x 10~ 5 eV 2 sin 2 #i2 — 0.3 [27], sin 2 #23 — 0.5,sin 2 #i3 — 
0.0 and Am| 3 = 2.1 x 10~ 3 eV 2 . Matter effects within the Earth are taken into account. Negative cosG^ corresponds to 
upward-going neutrinos and is the horizon. The electron neutrino flux ratio $ e /$e is shown in the right side of the figure. 
An expected excess (deficit) for atmospheric mixing in the first (second) octant is shown in the upper (lower) panel. The island 
shapes are regions of probability driven by the solar oscillation parameters. The center panel shows no significant region of 
excess or deficit when sin 2 #23 = 0.5. 



where P ex is the two neutrino transition probability 
(y e —> v x ) driven by Ara 2 2 and #12, L is the neutrino 
pathlength, and E is its energy. Using these equations 
the modified atmospheric v fluxes at Super-K become: 



$° [l + P ex (rcos 2 # 23 - 1)] 

COS 2 #23 



(5) 



(rCOS 2 #23 - 1) Pe: 



$0 



Sin 2 26>23 - V 1 ~ PexCOSCf 



where and <fr® are the neutrino fluxes in the absence 
of oscillations and r is their ratio. 

The left panel in Fig. [I] shows the transition probabil- 
ity, P exi as a function of energy and zenith angle, 9^, 
for neutrinos traversing the Earth (see below) assuming 
Am 2 2 = 7.7 x 10~ 5 eV 2 and sin 2 <9i 2 = 0.30 [27]. However, 
the overall effect of P ex on the electron neutrino flux at 
the detector is modified by the factor r cos 2 #23 — 1 as seen 
m Eq. (|J). Since the atmospheric neutrino flux ratio is 
r ~ 2 at low energies, there is no change in the v e flux 
if r cos 2 #23 = 1 (#23 = 7r/4). If cos 2 #23 is greater (less) 
than 0.5 (#23 < (>)tt/4) there is an expected enhance- 
ment (reduction) of the flux. Therefore it may be possible 
to determine the octant of #23 by observing changes in 
the flux of the low energy electron-like (e-like) samples 
at SK. Analogous changes to the flux on the other 
hand are suppressed by the leading factor of 1/r in the 
second term of Eq. Further, vacuum oscillations of 



the low energy flux are already well averaged so the 
correction appearing in the third term is negligible. The 
expected change in the v e flux as a function of energy 
and zenith angle for different values of #23 is shown as 
the right panel of Fig. [I] 

When #13 is different from zero, the matrix U13 is 
no longer sufficiently close to unity, and the above re- 
lations do not hold. Instead, in the search for non-zero 
#13, the oscillation analysis is done using a "one mass 
scale dominant" scheme wherein the solar neutrino mass 
difference is taken to be much smaller than the atmo- 
spheric mass difference. Accordingly, the solar mass dif- 
ference is neglected and a single mass splitting is adopted, 
Am 2 = m\ — m\ 2 such that Am 2 > 0(Am 2 < 0) cor- 
responds to the normal (inverted) mass hierarchy. In 
vacuum: 

P{y e v e ) 



1 - sin 2 2#i 3 sin 2 
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sin 2 #23sin 2 2#i3sin 2 
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Under this framework the three- neutrino oscillation prob- 
ability in constant density matter may be written [29] as 



P{yn v e ) = sin 2 



9 23 sin 2 2#^ 



1.27Am^L N 
E 
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FIG. 2: (color online). The three flavor oscillation probability Vy, -<-»• v e for #13 at the Chooz limit for 
neutrinos under the normal hierarchy in the one mass scale dominant framework is shown at left. In the 
right panel the v e flux ratio ^ hooz — l for oscillations with 613 at the Chooz limit relative to those at 
#13 = 0. Large matter-induced resonances between 2-10 GeV appear for upward going neutrinos traversing 
the core (cosG^ < —0.84) and mantle regions (—0.84 < cosG^ < —0.45). Atmospheric mixing is assumed at 
Amis = 2.1 x 10" 3 eV 2 and sin 2 6> 2 3 = 0.5. 



The matter modified mixing parameters are 



Am 2 



M 



2q/)M 



sin 2 2<9i 3 + (r-cos2<9i 3 ) 2 (8) 



sin 2 2(9 
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sin^2<9i3 + (r-cos2<9i 3 ) 2 



where T = ±.2y/2G fn e E / Am 2 , Gf is the Fermi con- 
stant, n e is the local electron density and the plus (mi- 
nus) sign specifies neutrinos (anti- neutrinos). Resonant 
enhancement of the oscillation probability occurs when 
|T| = cos 2#i3 and holds for either neutrinos or anti- 
neutrinos, depending on the mass hierarchy. Further, 
when #13 = there is no enhancement. 

Oscillation probabilities for neutrinos traversing the 
Earth appear in the left panel of Fig. [2J For Am 2 ~ 
2 x 10 -3 eV 2 this resonance occurs in the 2-10 GeV re- 
gion and its strength increases with #13 reaching ~ 40% 
conversion probability near the Chooz limit. Under these 
conditions, the primary signature in the atmospheric neu- 
trino sample at Super-K is an increased rate of high en- 
ergy upward-going e-like events. The right panel of the 
figure shows the v e flux ratio at SK oscillated with #13 at 
the Chooz limit relative to that oscillated at #13 = 0. Ad- 
ditional #i3-induced effects on muon event rates are ex- 
pected, but are generally much smaller. For large values 
of #13 an expected ~ 20% increase in the multi-ring e-like 
event (see below) rate would be accompanied by a ~ 5% 
change in similarly energetic muon-like (/x-like) samples. 

Including solar oscillation terms changes the oscillation 
probability in the resonance region by less than 5%, sup- 



porting our assumption of the "one mass scale dominant" 
framework. Their inclusion as an additional scanning pa- 
rameter also introduces a large computational burden in 
the #13 analysis and further motivates a separate #23 oc- 
tant analysis. 

Oscillation probabilities in both analyses are computed 
using a numerical technique [30] . Probabilities inside the 
Earth are computed using a piecewise constant radial 
matter density profile constructed as the median density 
in each of the dominant regions of the PREM [31 model: 
inner core (0 < r < 1220km) 13.0 g/cni , outer core 
(1220 < r < 3480km) 11.3 g/cm 3 , mantle (3480 < r < 
5701km) 5.0 g/cm 3 , and the crust (5701 < r < 6371km) 
3.3 g/cm 3 . Transition amplitudes are computed across 
each layer a neutrino traverses and the product of these 
together with the amplitude for crossing the Earth's at- 
mosphere is used to compute the final oscillation proba- 
bility. The difference in the obtained probabilities using 
this simplified model compared to the more expansive 
PREM model have a negligible impact on the final anal- 
ysis results after incorporating detector resolution effects. 



III. DATA SAMPLE 

Super-Kamiokande is a cylindrical 50 kton water 
Cherenkov detector situated at a depth of 2700 meters 
water equivalent. The detector volume is optically sep- 
arated into an inner volume (ID) and an outer veto re- 
gion (OD). During the SK-I (SK-II) periods the ID was 
instrumented with 11,146 (5,182) inward-facing 20-inch 
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FIG. 3: (color online). The distributions used in the tt° selection for five momentum regions: (a)P e < 250 MeV/c, 
(b)250 MeV/c < P e < 400 MeV/c, (c)400 MeV/c < P e < 630 MeV/c, (d)630 MeV/c < P e < 1000 MeV/c and 
(e)1000 MeV/c < P e . Solid (dashed) lines represent CCQE (NC) events in the FC sub-GeV single-ring e-like Monte 
Carlo. Events with an invariant mass above 100 MeV/c 2 are selected as 7r°-like. To separate 7r°-like and electron-like 
more efficiently, an additional likelihood selection is applied for events with momentum above 250 MeV/c. The 
distributions of the three likelihood variables are shown: the fraction of energy carried by the second fitted ring 
(E2KE1 + .E2)), the tt° mass and A-likelihood (described in the text). All distributions have been normalized to unit 
area. 



photomultiplier tubes (PMTs) and the OD with 1,885 
outward-facing 8-inch PMTs. In SK-III there were 11,129 
ID PMTs. Since SK-II, the ID PMTs have been encased 
in fiber-reinforced plastic shells with acrylic covers to pre- 
vent chain reactions within the detector in the event of 
a PMT implosion. A more detailed description of the 
detector may be found in [ 32] . 

In this paper, atmospheric neutrino events are orga- 
nized into three classes: fully contained (FC), partially 
contained (PC), and upward-going muons (UP/i). Events 
which deposit all of their Cherenkov light in the ID are 
classified as FC, while events that originate in the ID 
but have an exiting particle depositing energy in the OD 
are considered PC. Neutrino interactions occurring in the 
rock beneath the detector which produce muons that tra- 
verse the detector (through-going) or stop in the detector 
(stopping) are classified as UP/i events. Data accrued in 
the five years spanning the SK-I run period starting in 
1996 correspond to 1489 live-days of FC and PC events 
with 1646 UP/i live-days. SK-II data were taken be- 
tween December 2002 and October 2005 and represent 
799 (518) live-days of FC and PC events and 828 live- 
days of UP/i events. SK-III data were taken between 



December 2005 and June 2007 where the FC and PC 
livetime was 518 days and that for UP/i was 635. The 
difference of livetimes between FC/PC and UP/i is due 
to the insensitivity of the UP/i reduction to noise such as 
"flasher" PMTs. Such noise may be misconstrued as real 
FC/PC events so in those reductions data surrounding 
these events are rejected. 

Fully contained events are further divided into sub- 
GeV and multi-GeV sub-samples based on visible energy, 
E vis . Events with E vis < 1.33 GeV are considered sub- 
GeV. The number of reconstructed Cherenkov rings in an 
event is also used to separate these samples into single- 
and multi-ring sub-samples. Single-ring events are clas- 
sified into /i-like and e-like samples by the ring pattern. 
For multi-ring samples, the most energetic ring is used 
to classify the event type. Partially contained events are 
classified as "OD stopping" or "OD through-going" based 
on their energy deposition in the OD [33]. Similarly, 
UP/i events that traverse the detector are separated into 
"showering" and "non- showering" based on the method 
described in [34 while those that enter and stop within 
the detector are classified as "stopping." These samples 
are defined for all of the SK run periods. To enhance 
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FC sub-GeV single- 


-ring e-like 


FC sub-GeV 




0-decay 


1-decay 


7r°-like 


single-ring e-like 


MC events 


2663.2 


210.9 


191.8 


2996.4 


Q.E. 


77.7% 


3.8% 


10.6% 


70.6% 


CC single meson 


12.4% 


50.3% 


i .U /o 


10. Z /o 


v e + v e multi 7r 


1.0% 


9.7% 


1.8% 


1.7% 


coherent 7r 


1.3% 


8.5% 


0.5% 


1.7% 


CC Ufj, + 


0.6% 


15.2% 


7.0% 


2.0% 


NC 


6.8% 


11.2% 


72.0% 


8.7% 




FC sub-GeV single- 


■ring /x-like 


FC sub-GeV 




0- decay 


1-decay 


2-decay 


single-ring /x-like 


MC events 


1412.4 


2745.4 


164.3 


4297.8 


Q.E. 


71.3% 


78.5 % 


5.8% 


74.7% 


CC single meson 


12.9% 


15.5% 


65.7% 


16.7% 


Vfj, + Vfi multi 7r 


1.1% 


1.5% 


14.9% 


1.9% 


coherent tt 


0.8% 


1.5% 


8.6% 


1.6% 


CC zy e + u e 


1.8% 


< 0.1% 


< 0.1% 


0.7% 


NC 


11.8% 


2.6% 


3.3% 


4.3% 



TABLE I: The number of FC sub-GeV MC events and their fractional composition by neutrino interaction mode in SK-L The 
upper (lower) table shows the e-like (/i-like) sample. The left (right) side of the table shows the result after (before) separation 
into sub-samples. After separation, the CCQE purity is increased and the NC backgrounds are reduced in the 0-decay e-like 
and 1-decay /x-like sub-samples. 



each analysis' sensitivity to the desired oscillation effect, 
the FC samples have been further divided as outlined be- 
low. However, all of the data samples are used in both 
analyses. 

Several improvements to the reconstruction and Monte 
Carlo since earlier publications [lOj [TT] are incorporated 
in this paper. The ring counting likelihood has been up- 
dated to improve separation between the single-ring and 
multi-ring samples. Additionally, the neutrino interac- 
tion generator has been updated to include lepton mass 
effects in charged current (CC) interactions [35, 36 . An 
axial vector mass of 1.2 GeV has been used for quasi- 
elastic and single meson production processes and cross 
sections for deep inelastic scattering are computed based 
on the GRV98 parton distribution functions [37 . The at- 
mospheric neutrino flux is taken from [38] . More detailed 
information on the MC simulation, event generator, and 
event reconstruction is presented in [TT] . 



A. Additional sample selection for the 623 octant 
analysis 

To increase the purity of the interaction mode, FC sub- 
GeV single-ring events are separated into sub-samples 
based on their number of decay electrons and how 7r°-like 
they are. 

The FC sub- GeV single-ring e-like sample contains 
background events which are mainly neutral- current 
(NC) 7T° events where one of the two 7 rays from the 7r° 
decay has been missed by the event reconstruction. The 
electromagnetic shower from the 7 gives a light pattern 
similar to that of an electron and results in an electron- 



like classification. To reduce this type of background, a 
specialized 7r° fitter is used [39]. This fitter enforces a 
second ring on the data and then predicts a light pat- 
tern that would result from 7 rays propagating through 
the tank with the direction and vertex of the fitted rings. 
The intensity of each fitted ring as well as its direction 
are varied until the predicted light pattern best agrees 
with the observed one. Since the interaction mode of 
interest, charged-current quasi-elastic (CCQE), creates 
only one light-emitting particle, constructing the invari- 
ant mass for the two fitted rings provides some separation 
between CCQE and NC events. The left five panels in 
Fig. [3] show the invariant mass distributions from this 
tt° fitter for CCQE and NC events in the FC sub-GeV 
single-ring e-like Monte Carlo in five energy regions. Neu- 
tral current events tend to form a peak close to the 7r° 
mass whereas CCQE do not. For events with electron 
momentum below 250 MeV/c, a cut at 100 MeV/c 2 is 
used to create a 7r°-like sample. This cut, however, is not 
sufficient for higher electron momenta so an additional 
likelihood selection is used, incorporating three variables: 
the 7T° invariant mass distribution, the fraction of the 
event's reconstructed momentum carried by the second 
ring, and the difference of two likelihood variables which 
result from a 7r°-fit and electron- fit. The distribution of 
these variables is shown in Fig. [3J The 7r°-like selection 
likelihood functions are defined as, 

£ = ^io g (rf( a;i ))-iog(rf(x i )), (9) 

i=l 

where Tf fa)(T? fa)) represents the CCQE(NC) events' 
probability distribution function (PDF) for the i th vari- 
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able with observable X{. 

After separating the 7r°-like sample, the remaining 
e-like events are divided into two categories, 0-decay 
which has no decay electrons and 1-decay which has one 
or more decay electrons. Since v e CCQE events are not 
expected to produce decay electrons, there is a large frac- 
tion of CCQE interactions in the 0-decay sample. For the 
FC sub-GeV single-ring /i-like sample, there are three 
categories using the number of decay electrons: 0-decay, 
1-decay, and 2-decay, corresponding to the number of 
decay electrons reconstructed in the event. Since these 
CCQE events produce a muon they are expected to have 
at least one decay electron. Details of the event composi- 
tion by interaction mode after these event selections are 
shown in Tbl.|TJ The fraction of CCQE events is increased 
in the 0-decay e-like and 1-decay /i-like samples so they 
should improve sensitivity to changes in the sub-GeV flux 
induced by solar oscillations. 



B. Additional sample selection for the #13 analysis 

To improve sensitivity to v e appearance induced by 
non-zero #13, an enhanced FC multi-GeV multi-ring 
electron-like sample is created. The selection is based on 
a likelihood method [10] for SK-I and is extended in this 
analysis to SK-II and SK-III. The likelihood functions 
have been rebuilt using 100 years of MC incorporating 
recent improvements to the SK event reconstruction. Ac- 
cordingly, the event populations of the SK-I sample here 
differ from those in the reference. 

The MC is divided into five energy bins, 1.33-2.5 GeV, 
2.5-5 GeV, 5-10 GeV, 10-20 GeV, and > 20 GeV and 
PDFs for each bin are constructed using events whose 
most energetic ring has been reconstructed as electron- 
like. Four observables are used in the event selection: the 
number of decay electrons in the event, the maximum 
distance between the neutrino vertex and any muon de- 
cay electrons, the fraction of momentum carried by the 
event's most energetic ring, and the PID likelihood value 
of that ring. The final likelihood functions are defined 
as, 



where Yi represents the PDF for the i th observable and 
Xi is the observable's measured value. The superscripts 
S and B label the signal and background PDFs, respec- 
tively. The index j specifies the likelihood corresponding 
to one of the five energy bins considered. In selecting elec- 
tron neutrino events, the signal is taken to be CC u e -\-u e , 
while the background is composed of both CC v^+v^ and 
NC events. An event makes it into the final multi-GeV 
multi-ring sample if it passes all cuts in the FC reduc- 
tion, if the event's most energetic ring is electron-like, 
and if Cj > 0. Distributions of the likelihood variables 
for signal and background events appear in Fig. [4) De- 
cay electrons are produced in the signal sample through 
the decay chain of pions produced in these events. How- 
ever, lacking an exiting muon at the interaction vertex, 
fewer decay electrons are expected in the signal sample. 
Similarly, the maximum distance to a decay electron in 
the CC component of the background is expected to be 
larger due to the presence of energetic muons. The dis- 
tribution of the momentum fraction carried by the most 
energetic ring tends to peak towards higher values for 
signal events where the outgoing electron has been cor- 
rectly identified as electron-like. Background events, on 
the other hand, tend to peak at lower momentum frac- 
tions where the most energetic ring comes from a me- 
son or muon that has been misidentified as electron-like. 
Applying the likelihood improves the signal purity from 
53% to 74% in SK-I with 16% of the sample coming from 
NC events. Table [TT] shows the event compositions of 
the multi-GeV e-like sample after this selection for SK-I, 
SK-II, and SK-III. 

A summary of all atmospheric neutrino event samples 



used in this paper is shown in Tbl. Ill 



c d = ^iog(rf(x,))-iog(rf(x0), 



(10) 



IV. OSCILLATION ANALYSIS 

The oscillation analyses have been performed using the 
above data samples. Since the physical detector configu- 
ration differs between SK-I, SK-II, and SK-III, separate 
500 years-equivalent MC data sets for each run period 
are used. The data are compared against the MC expec- 
tation using a "pulled" \ 2 @Q] method based a Poisson 
probability distribution: 



2=1 



E n (l + J2f l n £i)-On + O n ln 



\ 



(11) 



where n indexes the data bins, E n is the MC expecta- bin. Systematic errors are incorporated into the fit via 
tion, and O n is the number of observed events in the n th the systematic error parameter where i is the system- 
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CC i/ ( 




CC ^ 




NC 


Total 




No £ 


c 


No £ 


C 


Nn C 

l\l U >C/ 


r 

JL, 




r 

JL, 


SK-I 


472.1 


331.0 


201.7 


39.2 


218.1 


74.4 


891.9 


444.6 


Percentage (%) 


53.0 


74.5 


22.6 


8.8 


24.5 


16.7 


100.0 


100.0 


SK-II 


253.4 


178.0 


110.7 


23.4 


119.1 


42.8 


483.2 


244.2 


Percentage (%) 


52.4 


72.9 


23.0 


9.6 


24.6 


17.5 


100.0 


100.0 


SK-III 


157.6 


112.8 


72.7 


15.1 


74.0 


26.1 


304.3 


154.0 


Percentage (%) 


51.8 


73.2 


23.9 


9.8 


24.3 


17.0 


100.0 


100.0 



TABLE II: The expected number of events for each interaction component of the multi-ring multi-GeV e-like sample before 
(No C) and after (£) likelihood selection for the SK-I, SK-II, and SK-III MC scaled to 1489, 798, and 518 livetime days, 
respectively. Two flavor neutrino oscillations ^ v T have been assumed with Am 2 = 2.1 x 10 ~ 3 eV 2 and sin 2 20 = 1.0. 





Data 


SK-I 

MC 


(osc.) 


Data 


SK-II 

MC 


(osc.) 


Data 


SK-III 

MC 


(osc.) 


FC Sub-GeV 




















single-ring 
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0- decay 


2984 


2655.9 


(2652.4) 


1605 


1405.8 


(1403.4) 


1098 


935.7 


(934.7) 


1-decay 


275 


204.4 


(194.3) 


155 


113.5 


(107.1) 


106 


69.6 


(66.7) 


7r°-like 


167 


159.1 


(155.2) 


81 


81.3 


(79.1) 


46 


45.6 


(44.4) 


//-like 




















0- decay 


1036 


1385.6 


(973.0) 


563 


765.9 


(537.2) 


346 


497.5 


(350.4) 


1-decay 


2035 


2760.6 


(1846.8) 


1043 


1429.4 


(957.3) 


759 


999.8 


(668.3) 


2- decay 


15U 


loo. 7 


(114.0) 




8z.8 


/r7 >-7\ 

(57.7) 


ol 


CO c 

58.5 


(41.0) 


2-ring 7r°-like 


497 


460.0 


(456.1) 


267 


237.3 


(235.1) 


178 


157.8 


(156.5) 


FC Multi-GeV 




















single-ring 




















e-like 


829 


777.8 


(777.7) 


392 


409.9 


(411.1) 


282 


279.3 


(278.4) 


//-like 


694 


1027.4 


(744.4) 


394 


550.2 


(399.0) 


231 


352.8 


(255.7) 


multi-ring 




















e-like 


433 


457.9 


(458.9) 


260 


252.3 


(251.9) 


149 


159.2 


(159.1) 


//-like 


617 


882.4 


(660.9) 


361 


459.6 


(344.1) 


226 


313.8 


(234.4) 


PC 




















OD stopping 


163 


222.7 


(167.3) 


116 


105.8 


(80.9) 


63 


75.1 


(55.7) 


OD through-going 


735 


965.4 


(755.0) 


314 


482.5 


(374.7) 


280 


334.9 


(262.8) 


Upward-going muon 




















stopping 


435.9 


701.7 


(419.4) 


207.6 


355.2 


(212.5) 


193.7 


273.8 


(163.5) 


non-showering 


1577.4 


1548.0 


(1343.9) 


725.3 


767.6 


(668.7) 


612.9 


599.4 


(520.0) 


showering 


271.6 


302.7 


(292.2) 


108.1 


147.8 


(143.6) 


110 


116.2 


(112.3) 


Reduction Efficiency 




















FC 






97.6% 






99.2% 






98.5% 


PC 






81.0% 






74.8% 






88.6% 


Upward stopping ji 






98.0% 






97.0% 






98.2% 


Upward through-going 


; m 




99.4% 






98.1% 






99.4% 



TABLE III: Summary of atmospheric neutrino data and MC event samples for FC, PC, and UP/x in SK-I, SK-II, and SK-III. 
The FC and PC livetime is 1489 days in SK-I, 798 days in SK-II, and 518 days in SK-III. The livetime of the UP/i samples is 
1645 days in SK-I, 827 days in SK-II, and 635 days in SK-III. The number of MC events has been normalized by the livetime of 
the data. The oscillated MC has been calculated using two flavor mixing at Am 2 = 2.1 x 10 -3 eV 2 and sin 2 20 = 1.0. 
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FIG. 4: (color online). Variables used in the likelihood definition to create the SK-I multi-GeV multi-ring e-like 
sample. The energy bins correspond to the the most energetic ring (MER) and the distributions have been scaled to 
the SK-I livetime. Signal events (CC v e + v e ) are shown as the solid line and background events (CC + Vy, and 
NC) are shown as the dashed line. The shapes of these distributions do not differ appreciably among SK-I, SK-II, 
and SK-III. 



atic error index and p n is the fractional change in the 
MC expectation in bin n for a 1-sigma change in the i th 
systematic error. The 1-sigma value of a systematic error 
is labeled as c^. Equation (11) is minimized with respect 



to the €i at each point in a fit's oscillation parameter 
space according to = 0. This derivative yields a set 
of linear equations in e that can be solved iteratively [40] . 
The best fit point is defined as the global minimum y 2 
on the grid of oscillation points. 



To ensure stability of the function in Eq. (11) the bin- 



ning has been chosen so that there are at least 6 expected 
MC events in each bin after scaling to the SK-I livetime. 
Data are binned separately for SK-I, SK-II, and SK-III, 
each with a total of 420 bins. Both analyses simultane- 
ously fit 16 event samples, including both e-like and /i-like 
event categories, shown in Tbl. |III| The samples sepa- 
rated by number of decay electrons are divided into 5 mo- 
mentum and 10 zenith angle bins for the 0-decay e-like, 
0-decay and 1-decay /i-like samples, and 1 zenith bin oth- 
erwise. The remaining FC events are divided among 14 
momentum bins, PC events into a total of 6 bins, and all 
upward through-going muon samples have one momen- 
tum bin each. The upward-stopping muon samples have 
been divided into three momentum bins. All of these 
samples are further divided into 10 evenly-spaced zenith 



angle bins. FC and PC events range from — 1 < cos0 < 1 
and UP/i events are binned from — 1 < cosO < 0. 

Both analyses consider 120 sources of systematic un- 
certainty. These systematic errors are separated into two 
categories, those that are common to all of the SK run 
periods and those that differ. Errors that are classi- 
fied as common are related to uncertainties in the atmo- 
spheric neutrino flux, neutrino interaction cross sections, 
and particle production within nuclei. Systematic errors 
that are independent for SK-I, SK-II, and SK-III repre- 
sent uncertainties related to the detector performance in 
each era. Particle reconstruction and identification un- 
certainties, as well as energy scale and fiducial volume 
uncertainties, differ for SK-I, SK-II, and SK-III because 
of their different geometries. These systematics are there- 
fore considered as separate sources of uncertainty. The 
effect of the systematic uncertainties are introduced by 
the coefficients f % n which are computed for every bin and 
error in the analysis. For common systematic uncertain- 
ties there is a coefficient for every bin in the analysis. On 
the other hand, independent systematic errors specific to 
SK-I (II, III) have non-zero coefficients for the SK-I (II, 
III) analysis bins and are zero otherwise. Tables [V] and 
[Vljlist the 33 common errors separated into neutrino flux 
and interaction-related systematics, respectively. Table 



VII lists the 29 x 3 independent systematic errors and all 
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three tables show errors with their fitted value, e^, from 
the #13 search, together with their uncertainty. More in- 
formation about these systematic errors is presented in a 
previous analysis [TT] . 

To prevent instabilities in the y 2 value resulting from 
the low statistics data in later SK run periods, the SK-II 
and SK-III bins are merged with those of SK-I. In the 
minimization of the function in Eq. (11) the following 
changes are made: 

O n -»■ J2°n Kl 

i 

E n (i+J2ffc) -»• £i£ m (i +£/&■)■ 




0.35 



0.4 0.45 0.5 0.55 0.6 0.65 



Since the systematic error coefficients are computed for 
separate SK-I, SK-II, and SK-III bins as discussed above, 
merging in this way preserves the effect of the system- 
atic errors specific to each detector geometry. Using this 
method, the final y 2 is taken over 420 merged bins. 



FIG. 5: (color online). X 2 ~Xmin distribution as a function 
of sin 2 O23 for oscillations without the 1-2 parameters (dotted 
line) and with the 1-2 parameters (solid line). For each sin 2 #23 
point, Am23 is chosen so that x 2 is minimized. The horizontal 
line corresponds to the 68% (90%) confidence level which is 
located at xLin + 1-0(2.7). 



A. 623 octant analysis 

In the search for 623 ^ 7r/4, two fits are performed 
to the data to extract a constraint on sin 2 O23 assum- 
ing #13 = 0. The first (solar-off) is done over the 
two-dimensional space of Am^ and sin 2 # 23 (41 points 
each of Ara| 3 in [1.0,6.3] x 10~ 3 eV 2 and sin 2 <9 23 in 
[0.3,0.7]) and is compared to a second (solar-on) fit, ex- 
panded to four dimensions including the solar oscilla- 
tion parameters Ara 2 2 and sin 2 #i 2 (fit over 4 points of 
Aml 2 in [7.41,7.94] x 10" 5 eV 2 and 5 points of sin 2 #i 2 in 
[0.28,0.36]). This grid of points has been chosen based 
on a combined fit of the solar neutrino experiment and 
KamLAND data [271 SB- To constrain the fit over the 
solar parameters, the ^xloiar va l ue from this combined 
analysis is then added to that of the fit at each of these 
grid points. 

Figure [5] shows the A% 2 distributions with and with- 
out the solar parameters as a function of sin 2 # 2 3, where 
Ara 2 2 , Ara 23 , and sin 2 #i 2 are chosen so that y 1 is mini- 
mized. The best-fit point with the solar parameters is 



located at sin O23 — 0.50, 



Aml 3 



2.1 x 10~ 6 eV z 



12 = 0.30, and Amf 2 = 7.59 x 10" 5 eV z , with a min- 
2 of 470.6/416 d.o.f., while that without the solar 

2tf _ n A™2 o 1 w in— 3 n \r2 



imum x 

parameters is sin 2 ^ = 0.50, Am 23 = 2.1 x 10 -3 eV 
with a minimum x 2 of 469.6/418 d.o.f. No significant de- 
viation of sin 2 # 2 3 from 7r/4 is seen with the addition of so- 
lar terms to the analysis but they do give rise to the asym- 
metric shape seen in the x 2 distribution. Including the 
solar terms constrains the measurement of sin 2 # 2 3 at the 
68 (90)% C.L. to 0.438(0.407) < sin 2 <9 23 < 0.558(0.583). 
The up-down asymmetry of the single-ring e-like data in 
comparison with the best fit MC expectation and the ex- 
pectations for sin 2 # 2 3 = 0.4 and 0.6 appears in Fig [6] 
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FIG. 6: (color online). Asymmetry (Up - Down) /(Up + 
Down) for the SK-I+II+III single-ring e-like 0-decay data set 
in the octant analysis. Up is defined as events with cosB < 
—0.2 and down as cosG > 0.2. The solid line represents the 
MC expectation at the best fit point and the dashed (dotted) 
line shows the expected asymmetry for sin 2 #23 = 0.4(0.6). 
The error bars are statistical. 



B. 0i3 analysis 

In the #13 analysis, oscillation fits are performed by 
scanning a grid of 83,025 oscillation points in three 
variables: log 10 Am 2 , sin 2 023? and sin 2 #i3. The fitting 
procedure has been performed on the combined SK- 
I+II+III data set assuming both a normal and inverted 
hierarchy. The best fit in the normal hierarchy is at 
Am 2 = 2.1 x 10- 3 eV 2 , sin 2 <9i 3 = 0.0, and sin 2 <9 23 = 0.5 



with x 



mm 



468.7. In the inverted hierarchy, the fit 
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SK-I+II+III sin 2 / 


9l3 


90% C.L. 


sin 2 #i 3 


Am 2 


[eV 2 ] 


sin 2 #23 


x 2 /do.f 


Normal Hierarchy 


< 


0.04 


0.00 


2.1 x 


10" 3 


0.50 


468.7 / 417 


Inverted Hierarchy 


< 


0.09 


0.006 


2.1 x 


10" 3 


0.53 


468.4 / 417 



TABLE IV: Best fit information for fits to the SK-I+II+III data for both hierarchies in the 613 analysis. The limit on sin 2 613 
is the C.L. in one-dimension at 90% and the corresponding bounds on Am 2 are 1.9(1.7) x 10~ 3 < Am 2 < 2.6(2.7) x 10~ 3 eV 2 
in the normal (inverted) hierarchy. 
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FIG. 7: (color online). Normal hierarchy allowed regions at 68% (thin line), 90% (medium ), and 99% (thick) C.L. 
for the SK-I+II+III data. The shaded region in the first panel shows the Chooz 90% exclusion region. 



is at 
sin 2 (9. 



23 



Am 2 = 2.1 x 10~ 3 eV^ s'm 2 6 13 = 0.006, and 
= 0.53. The results are summarized in Tbl. IV 



No preference is seen in the data for either mass hier- 
archy. Confidence intervals at 90% (99%) are drawn in 
two dimensions at y 2 = y 2 • + 4.6(9.2) and the third 
parameter point in these projections has been minimized 
over at each point in the plane to give the smallest value 
of x 2 ' Computing the 90% (99%) critical value using a 
Feldman-Cousins [42] type procedure confirmed that 4.6 
(9.2) is the correct value in this scheme. The resulting 
allowed regions and corresponding A% 2 distributions are 
shown in Figs.[7| [8j and[9j In the first panel of the former 
two figures, the Am 2 vs. 6 is plane for the fit is overlaid 
with the Chooz [23] 90% C.L. exclusion contour. The 
zenith angle distributions of the combined data overlaid 
with the best fit MC expectation in the normal hierar- 
chy and the expectation resulting from #13 at the Chooz 



limit are shown in Fig.[T0j Figure [TT] shows the up-down 
asymmetry for the single-ring (multi-ring) e-like sample 
as a function of lepton momentum (total energy) for the 
single-ring (multi-ring) sample. The data are consistent 
with sin 2 #i3 = 0. 



V. CONCLUSION 

A three flavor oscillation fit to the first, second, and 
third generation Super-K atmospheric neutrino data has 
been performed. No evidence for #13 > is found in 
fits to either hierarchy. The best fit oscillation parame- 
ters in the normal (inverted) hierarchy are Am 2 = 2.1 x 
10- 3 eV 2 , sin 2 0i3 = 0.0(0.006), and sin 2 <9 23 = 0.5(0.53). 
The value of 6 is is constrained to sin 2 #13 < 0.04(0.09) 
at the 90% confidence level. All fits are consistent with 
the Chooz experiment's upper limit and no preference for 
either mass hierarchy exists in the data. The #23 octant 
analysis finds no evidence for a preferred octant for #23- 
However, the mixing angle is constrained at 90% C.L. to 
0.407 < sin 2 <9 23 < 0.583. 
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FIG. 9: (color online). The A% 2 distributions for the SK-I+II+III data set in the normal (solid) and inverted 
(dashed) hierarchies. The horizontal lines represent the position of the one-dimensional cut value corresponding to 
the (from top to bottom) 99%, 90%, and 68% CL. 
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FIG. 10: (color online). SK-I+II+III zenith angle and lepton momentum distributions of the event samples used in the 
analyses. Black dots represent the data with statistical errors, the solid lines are the MC expectation at the best fit from the 
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A. Appendix 



Tables [Vj [VlJ and |VII| summarize the best fit system- 
atic error parameters for the best fit point in the normal 
hierarchy fit from the #13 analysis. 
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Systematic Error 






fit value 


a 


Flux normalization 


E u < 1 GeV 




34.7 25 a 




E u > 1 GeV 




8.8 


7 b 


H" 1 c 












E u < 1 GeV 




-1.9 


2 




1< E v < 10 GeV 




-2.5 


3 




E u > 10 GeV 




-3.7 


5 C 














E u < 1 GeV 




5.54 


5 




1< Eu < 10 GeV 




1.13 


5 




> 10 GeV 




-0.10 


8 d 














Eu < 1 GeV 




-0.48 


2 




1< Eu < 10 GeV 




-1.35 


6 




> 10 GeV 




-1.75 


6 e 


Up /down ratio 


< 400 MeV 


e-like 


-0.07 


0.1 






/i-like 


-0.23 


0.3 






0- decay u-like 

J r 


-0.84 


1.1 




> 400 MeV 


e-like 


-0.61 


0.8 






/i-like 


-0.38 


0.5 






0- decay /x-like 


-1.29 


1.7 




Multi-GeV 


e-like 


-0.53 


0.7 






/i-like 


-0.15 


0.2 




Multi-ring Sub-GeV 


/i-like 


-0.15 


0.2 




Multi-ring Multi-GeV 


e-like 


-0.23 


0.3 






/i-like 


-0.15 


0.2 




PC 




-0.15 


0.2 


Horizontal/ Vertical ratio < 400 MeV 


e-like 


-0.01 


0.1 






/i-like 


-0.01 


0.1 






0- decay /i-like 


-0.03 


0.3 




> 400 MeV 


e-like 


-0.14 


1.4 






/i-like 


-0.19 


1.9 






0- decay /x-like 


-0.14 


1.4 




Multi-GeV 


e-like 


-0.33 


3.2 






/i-like 


-0.23 


2.3 




Multi-ring Sub-GeV 


/i-like 


-0.13 


1.3 




Multi-ring Multi-GeV 


e-like 


-0.29 


2.8 






/i-like 


-0.15 


1.5 




PC 




-0.17 


1.7 


K/tt ratio in flux calculation 




-12.9 


10' 


Neutrino path length 






-8.8 


10 


Sample- by-sample 


FC Multi-GeV 




-4.5 


5 




PC + Up-Stop jJL 




-7.1 


5 



Uncertainty linearly decreases with log^ from 25%(0.1GeV) to 7%(lGeV). 

^Uncertainty is 7% up to 10 GeV, linearly increases with log^ from 7%(10GeV) to 12%(100GeV) and then to 20%(lTeV) 
Uncertainty linearly increases with log E u from 5 %(30 GeV) to 30 %(1 TeV). 
^Uncertainty linearly increases with log^ from 8%(100GeV) to 20%(lTeV). 
Uncertainty linearly increases with log^ from 6%(50GeV) to 40%(lTeV). 
'Uncertainty increases linearly from 5% to 20% between lOOGeV and ITeV. 



TABLE V: Flux-related systematic errors that are common to all SK geometries. The second column shows the best fit value 
of the systematic error parameter a in percent and the third column shows the estimated 1-a error size in percent. 
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Systematic Error 




fit value 


a 


MA in QE and single tt 




-2.4 


10 


CCQE cross section 




0.66 


1.0° 


Single meson production cross section 




7.8 


20 


DIS cross section (E nu < 10 GeV) 




-0.16 


1.0 5 


DIS cross section 




2.27 


5 


Coherent tt production 




1.53 


100 


NC/(CC) 




1.51 


20 


Nuclear effect in 16 nucleus 




-13.8 


30 


Nuclear effect in pion spectrum 




0.8 


1.0 C 


v T contamination 




1.0 


30 


NC in FC /i-like (hadron simulation) 




-4.6 


10 


CCQE Vijvi (i=e,/x) ratio 




0.84 


1.0 a 


CCQE fi/e ratio 




1.12 


1.0 a 


Single tt production, tt° /7r ± ratio 




-29.0 


40 


Single tt production, Vijvi (i=e,/i) ratio 


-0.04 


1.0 d 


7r + decay uncertainty Sub-GeV 1-ring 


e-like 0- decay 


-0.48 


0.5 




/x-like 1- decay 


0.77 


-0.8 




e-like 1-decay 


3.9 


-4.1 




/x-like 0-decay 


-0.77 


0.8 




/x-like 2-decay 


5.46 


-5.7 



a Difference from the Nieves |43| model is set to 1.0 
^Difference from CKMT 44 parametrization is set to 1.0 
difference between NEUT [Jg] and NUANCE 46 is set to 1.0 
^Difference from the Hernandez |47| model is set to 1.0 

TABLE VI: Neutrino interaction and particle production systematic errors that are common to all SK geometries. The second 
column shows the best fit value of the systematic error parameter 6i in percent and the third column shows the estimated 1-a 
error size in percent. 
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SK-I SK-II SK-III 



Systematic Error 




fit 


a 


fit 


a 


fit 


a 


FC reduction 




0.005 


0.2 


0.008 


0.2 


0.061 


0.8 


PC reduction 




-0.99 


2.4 


-2.12 


4.8 


0.034 


0.5 


FC/PC separation 




-0.058 


0.6 


0.068 


0.5 


-0.28 


0.9 


PC-stop/PC-through separation (top) 

yr i o ± \ it / 




7.84 


14 


-17.47 


21 


-20.03 


31 


PC-stop/PC-through separation (barrel) 




-2.27 


7.5 


-31.51 


17 


3.44 


23 


PC-stop/PC-through separation (bottom) 




-2.32 


11. 


-7.32 


12 


1.59 


11 


Non-zy BG (e-like) Sub-GeV 




0.077 


0.5 


0.004 


0.2 


0.003 


0.1 


Multi-GeV 




0.047 


0.3 


0.005 


0.3 


0.011 


0.4 


Non-zy BG (/x-like) Sub-GeV 




-0.01 


0.1 


0.02 


0.1 


0.052 


0.1 


Multi-GeV 




-0.01 


0.1 


0.02 


0.1 


0.11 


0.2 


Sub-GeV 1-ring 


/i-like 0-decay 


-0.04 


0.4 


0.02 


0.1 


0.052 


0.1 


PC 




-0.02 


0.2 


0.14 


0.7 


0.95 


1.8 


Fiducial volume 




-0.23 


2 


0.43 


2 


0.93 


2 


Ring separation < 400 MeV 


e-like 


1.23 


2.3 


-1.67 


1.3 


0.12 


2.3 




/i-like 


0.37 


0.7 


-2.96 


2.3 


0.037 


0.7 


> 400 MeV 


e-like 


0.21 


0.4 


-2.19 


1.7 


0.021 


0.4 




/i-like 


0.37 


0.7 


-0.90 


0.7 


0.036 


0.7 


Multi-GeV 


e-like 


1.97 


3.7 


-3.35 


2.6 


0.19 


3.7 




/i-like 


0.91 


1.7 


-2.19 


1.7 


0.089 


1.7 


Multi-ring sub-GeV 


/i-like 


-2.40 


-4.5 


10.56 


-8.2 


-0.24 


-4.5 


Multi-ring multi-GeV 


e-like 


0.05 


0.1 


-2.45 


1.9 


0.16 


3.1 




/i-like 


-2.19 


-4.1 


1.03 


-0.8 


-0.21 


-4.1 


Particle identification Sub-GeV 


e-like 


-0.007 


0.1 


0.13 


0.5 


0.004 


0.1 




/i-like 


0.007 


-0.1 


-0.13 


-0.5 


-0.004 


-0.1 


Multi-GeV 


e-like 


-0.014 


0.2 


0.023 


0.1 


0.008 


0.2 




/i-like 


0.014 


-0.2 


-0.023 


-0.1 


-0.008 


-0.2 


Particle identification (multi-ring) Sub-GeV 


/i-like 


-0.18 


-3.9 


-0.55 


-2.2 


-0.15 


-3.9 


Multi-GeV 


e-like 


0.078 


1.7 


0.45 


1.8 


0.063 


1.7 




/i-like 


-0.13 


-2.9 


-0.86 


-3.4 


-0.11 


-2.9 


Energy calibration 




-0.002 


1.1 


-0.56 


1.7 


-0.35 


2.7 


Up/Down asymmetry energy calibration 




-0.4 


0.6 


-0.15 


0.6 


-0.03 


1.3 


Upward-going muon reduction Stopping 




-0.057 


0.7 


-0.14 


0.7 


0.14 


0.7 


Through-going 




-0.041 


0.5 


-0.10 


0.5 


0.10 


0.5 


Upward stopping/through-going \i separation 




-0.04 


0.4 


0.006 


0.4 


0.04 


0.6 


Energy cut for upward stopping \i 




-0.13 


0.8 


-0.26 


1.4 


0.78 


2.1 


Path length cut for upward through-going a 

O IT o O O / 




0.39 


1.8 


-1.0 


2.1 


0.4 


1.6 


Upward through-going a showering separation 




9.42 


9.0 


2.28 


13.0 


6.1 


6.0 


BG subtraction of upward /x a Stopping 




4.16 


16 


-7.47 


21 


0.004 


20 


Non- showering 




-1.24 


11 


8.08 


15 


6.34 


19 


Showering 




2.27 


18 


-18.16 


14 


24.7 


24 


Multi-GeV Single-Ring Electron BG 




5.95 


16.3 


-4.67 23.4 


1.06 41.4 


Multi-GeV Multi-Ring Electron BG 




-4.38 35.6 


-1.4 22.3 


-16.8 


38.0 


Multi-GeV Multi-Ring e-like likelihood 




-1.12 


6.4 


0.5 


11.1 


-0.3 


5.3 


Sub-GeV 1-ring tt° selection 100 < P e < 250 


MeV/c 


-3.94 


11.2 


-4.08 


7.5 


-5.34 


7.7 


250 < P e < 400 




-4.05 


11.5 


-4.85 


8.9 


-18.37 26.4 


400 < P e < 630 




-8.23 23.4 


-9.52 


17.5 


-8.70 


12.5 


630 < P e < 1000 




-6.72 


19.1 


-5.81 


10.7 


-18.58 26.7 


1000 < P e < 1330 




-4.57 


13.0 


-6.03 


11.1 


-18.58 26.7 


Sub-GeV 2-ring tt° 




-0.31 


2 


0.024 


2 


0.009 


1 


Decay-e tagging 




0.16 


1.5 


0.41 


1.5 


1.06 


1.5 


Solar Activity 




0.6 


20 


27.9 


50 


3.78 


20 



a The uncertainties in BG subtraction for upward-going muons are only for the most horizontal bins —0.1 < cos# < 0. 

TABLE VII: Systematic errors that are independent between SK-I, SK-II, and SK-III. Columns labeled "fit" show the best fit 
value of the systematic error parameter a in percent. Those labeled "a" show the estimated 1-a error size in percent. 



